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ABSTRACT 

In general, yaw control of aircraft plays a major role. It is controlled by rudder deflection by creating side force 
with help of rudder pedals. The aircraft yaw is controlled by several methods such as Fuzzy logic, linear quadratic 
controller (LQR), optimal pole placement method. These methodologies achieved only the transient response, but failed to 
obtain the steady state response. There are some cases, where the steady state response is necessary. In this research paper, 
we propose a new methodology called Linear Algebraic approach, where both the transient and steady state responses are 
achieved and the disturbance conditions are totally suppressed. 
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INTRODUCTION 

After the first flights, articulated control surfaces were introduced for basic control which is operated by the pilot 
through a system of cables and pulleys. This technique survived for decades and at present used for small airplanes. 
The solution to this problem was made by the introduction of aero dynamic balances and tabs. But further growth of the 
aircraft and flight envelopes enhanced the need of powered systems to control the articulated aerodynamic surfaces. 

Modern aircraft include a variety of automatic control system that aids the flight crew in navigation, flight 
management and augmenting the stability characteristic of the airplane. Designing an autopilot requires control system 
theory background and knowledge of stability derivatives at different altitudes and Mach numbers for a given airplane. 
Aircraft have a number of different control surfaces: the primary flight controls are pitch, roll and yaw controls which are 
controlled by deflection of elevators, ailerons and rudder. Yaw is controlled by the rudder. The pilot moves the rudder 
sideways and the necessary yaw angle to change the direction of the flight as per requirement of mission profile of the 
aircraft. 

A modern controller (LQR) and intelligent fuzzy logic controller (FLC) is Developed for control the yaw of an 
aircraft system, but it failed to suppress the disturbances. 

In this research paper, the control system design for yaw control is designed. A feed forward and feedback 
controller design is developed for control of the yaw of an aircraft system. In this the disturbances are totally suppressed 
based on the wind axis direction and flight path envelope. 

MODELING OF YAW CONTROL 

System Two types of dynamical equations are present for an aircraft. The lateral dynamic equations of motion, 
which represents the dynamics of aircraft with respect to lateral axis and longitudinal dynamic equations of motion which 
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represents the aircraft's dynamics with respect to longitudinal axis. Lateral dynamics includes yaw, roll and sideslip 
motions of aircraft. Pitching motion comes under longitudinal dynamics. In this research paper, control of yaw angle of 
aircraft, when it performs yawing motion is explained. The control surfaces of the aircraft are shown in Figure 1 . 




Figure 1: Yaw, Roll & Pitch Motion of Aircraft [1] 

The forces, moments and Velocity components in the body Fixed frame of an aircraft system are as shown in 
figure 2, where L, M and N represent the aerodynamic moment components, the term p, q and r represent the angular rates 
components of roll, pitch and yaw axis and the term u, v and w represent the velocity components of roll, pitch and yaw axis. 




Figure 2: Body Fixed Frame of Aircraft 

For deriving the lateral equations we assumed that the aircraft is in steady cruise with constant altitude and velocity. 
Also it is assumed that change in pitch angle does not change the speed of the aircraft and the reference flight conditions are 
symmetric with propulsive forces constant. Therefore, 



v=p=q=r= ip = 
Y+mgC e S e = m(^ + ru-pw) 



N = -hz^+Iz^ + P?Uy ~ *x} ~ hzV 



(1) 

(2) 

(3) 
(4) 



The equations are linearized using small- disturbance theory by replacing all the variables in equations (2), (3) and 
(4) with a reference value plus a small disturbance equation (5) 

u = u + Au; v = v + Av; w = vt T + Aw 

p = p +Ap;q = q +Aq;Y = Y + AY 

r = r Q + Ar;L = L + AL;M = M + AM 
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S = S a + AS 



(5) 



(6) 



(7) 



(8) 



The lateral directional equations of motion consist of the side force, rolling moment and yawing moment 
equations of motion (6), (7) and (8) represents the linearized form. We are taking sideslip angle Ap instead of the side 
velocity Av. These two quantities are related to each other in the following way, [3] 



A/? * tan" 1 « ^ 

Lateral equations of motion in state space form is shown in equation (10) 
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(10) 



For this system, the input will be the aileron deflection angle and the output will be the yaw angle. In this study, 
the data from NAVION Transport [2] is used in system analysis and modeling. The lateral directional derivatives stability 
parameters for this airplane are given in Table 1. The values are taken directly from reference [3], as these are standard 
data of the NAVION aircraft. 

Table 1: The Lateral Directional Derivatives Stability Parameters [3] 



General 
A^a ti <:» n 
Aii-plsa no: 
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Y-Force 
Deri\a t i\"e & 


Yawing 
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The values in the above table is substituted in equation (10). The aileron deflection S c in (10) is neglected as we 
are only concerned about rudder deflection 5 r . 
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(11) 



Transfer function from rudder deflection angle to yaw angle is given by the equation (12) 

A _ -4.6130s 3 - 47.95625= - ll.BS33s -f 5.7410 
A£-(s) ~ s 4 -h 9.4090s- + 14.01S9s" -f 43.4991s -f 0.3979 

3. Two -Parameter Controller Design 



(12) 









CDmpwHloi 


m 


CLE 



fatftMfc 



Figure 3: General Two Parameter Controller 

The general two-parameter controller configuration is shown in Figure 3. The two-parameter controller 
configuration used in this work is shown in Figure 3. Here G(s) is plant transfer function, Ci(s) and C 2 (s) are the two 
compensators, the two compensators will be chosen to have same denominator and have the form 



Ci(s) =L(s)/A(s). 
C 2 (s) =M(s)/A(s). 

Where L(s), M(s) and A(s) are polynomials to be determined. The transfer function from R to Y is 
G (s)= N(s) L(s)/ [A(s) D(s) +M(s) N(s)] 
= No(s)/D (s) 



(1) 



(2) 



In equation (2) if all the roots of D (s) can be arbitrarily assigned, then the design is said to achieve arbitrary pole 
assignment. Assume G(s) as N(s)/D(s) which is strictly proper. Also N(s) and D(s) are co-prime. Then as per [5] arbitrary 
pole assignment is always possible. If the degree of the compensator C 2 (s) < (n-1), arbitrary pole placement is not possible 
for every set of poles. The design procedure for the implementation of the overall transfer function Go (s) is as follows: 



Step 1: Compute the following rational function 
G (s) /N(s) = N (s)/{ D (s) N(s)}=tf B (s)/D p (s) 



(3) 



Impact Factor (JCC): 5.9638 



Index Copernicus Value (ICV):3.0 



Robust Controller for Aircraft Yaw 



69 



Where iV a (s) and I? a (s) are co-prime. If No(s) and Do(s) are co-prime then common factors may exist only 
between No(s) and N(s). Cancel all common factors between them by computing G (s)/N(s). 

Step 2: From equation (3) the degree of A(s) D(s) + M(s) N(s) is deg {A(s)D(s)+M(s)N(s)}= 
deg{A(s)}+deg{D(s)}=2n-l Since C L (s) and C = (s) are proper, deg N(s) < deg D(s) and deg M(s) < deg A(s). If deg 
Dp(s) = p < 2n-l, in order to get C 2 (s) as proper and to match the coefficients of the same power of V on both sides of the 
equation, introduce an arbitrary polynomial D a (s) of degree 2n-l-p to get the total degree Z? a (s) = (2n-l) = deg {A(s) D(s) + 
M(s) N(s)}. Because this polynomial will be cancelled in the design procedure, its roots should be chosen in the acceptable 
pole zero cancellation regions. The cancelled roots i.e., the roots of DL(s) must be chosen three to four times faster than the 
poles of Go(s).If the deg D a (s) =p = (2n-l), set D (s) =1. In most applications we have deg D a (s) < 2n-l, the case in which 
deg D (s) > 2n-l is not considered. If G(s) is bipolar i.e., deg N(s) equals to deg D(s), then the above procedure in step 2 is 
modified as if deg D a (s) = p < 2n, introduce an arbitrary polynomial E(s) of degree (2n-p). 

Step 3: G (s) = {N(s) JV p (s)}/ D a (s) 

=N(s){E(s)JV p (s)}/E(s)i^(s) (4) 
Comparing (2) and (4) we get 

L(s) = JV p (s)E(s) (5) 

A(s) D(s) + M(s) N(s) = E(s) D a (s) = F(s) (6) 

Equation (6) is called the Diophantine equation. In this equation D(s) and N(s) are known, the roots of D a (s) are 
known as closed loop poles to be assigned from N(s) and the roots of E(s) are arbitrarily chosen. The polynomials A(s) and 
M(s) are found. We solve (6) by matching the coefficients of the same powers of Von both sides of the equation i.e., we 
solve a set of linear algebraic equations as in [4] to get A(s), M(s) and L(s). Knowing these polynomials we can implement 
Go(s), C L (s) and C 2 (s). The introduction of E(s) in the design procedure is crucial. If we introduce E(s) as suggested in 
step 2, then M (s)/A(s) and L(s)/A(s) will be proper, and the resulting system is well posed. The design involves pole zero 
cancellations. The canceled poles are the roots of E(s), which are chosen by the designer. Thus, if G (s) is stable and E(s) is 
Hurwitz, then the system is totally stable. 

SIMULATION AND RESULTS 

Simulations are first carried out on unity feedback system with and without disturbance. It is observed from 
Figure 4 although the steady state error is zero for a step input, and %overshoot is Zero for compensated system. 

From Figure 5 we observe the compensated system totally suppressed disturbance at time 5 sec where as the 
Uncompensated system have large %over shoot and infinite steady state error 
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Figure 4: Response for Without Disturbance Condition 
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Figure 5: Response for With Disturbance Condition 

CONCLUSIONS 

Thus, studies are carried out on aircraft yaw control, by incorporating two parameter controllers with linear 
algebraic method. The total system rejects the disturbance by tracking the reference input with zero steady state error. 
In general the procedure yields better results. 
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